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PhI Abstract 

D , We adapt a general method to solve both the full and reduced Salpeter 

equations and systematically explore the conditions under which these two 
^ ■ equations give equivalent results in meson dynamics. The effects of constituent 

\^ • mass, angular momentum state, type of interaction, and the nature of confine- 

ment are all considered in an effort to clearly delineate the range of validity 
of the reduced Salpeter approximations. We find that for J 7^0 the solutions 
are strikingly similar for all constituent masses. For zero angular momentum 
states the full and reduced Salpeter equations give different results for small 
quark mass especially with a large additive constant coordinate space poten- 
tial. We also show that ^ corrections to heavy-light energy levels can be 
accurately computed with the reduced equation. 



1 Introduction 



The instantaneous Bethe-Salpeter equation, or Salpeter equation is by far the 
most commonly employed relativistic wave equation in meson models with fermionic 
constituents. Until recently, almost all explicit calculations have used a simplified 
version known as the reduced Salpeter equation. The later becomes identical to the 
full Salpeter equation if at least one of the constituent masses is infinite. 

The reduced Salpeter equation is of the standard eigenvalue (hermitian) type 
whereas the full equation is not. Its solutions are thus algebraically and numerically 
simpler than that of the full equation. For example, the reduced equation doesn't 
have negative energy solutions, nor does it have solutions with zero norm, both of 
which exist for the full Salpeter equation [0, ||. More importantly, the reduced 
equation has variationally stable solutions for a wider range of kernel types than 
does the full equation |^. For example, there are no variationally stable solutions 
to the full Salpeter equation corresponding to pure scalar confinement. The reduced 
Salpeter equation, on the other hand, has well defined variationally stable solutions 
with scalar confinement. Also, the reduced equation is equivalent to the "no-pair" 
equation 0] proposed to cure the "continuum dissociation" problem in relativistic 
atomic physics. There are therefore historical, practical, and physical reasons for 
using the reduced equation. We outline here the conditions under which this can be 
done without sacrificing accuracy. 

In the real world the constituent mass is never infinite, so one faces a quantitative 
question as to the practical region of validity of the reduced Salpeter equation. Our 
results here establish that for many purposes the reduced Salpeter equation is quite 
adequate and one can take advantage. An analysis involving heavy-light mesons with 
c or & quarks, or bb, cc, or ss onia states, does not incur serious error by using the 
reduced Salpeter equation. It is only for J = states and with small quark masses 
where there can be significant differences between the full and reduced Salpeter 
solutions. Dynamical models involving light pseudo-scalar states such as the vr, r], 
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or K mesons can lead to serious errors if the full Salpeter equation is not used. 

Our analysis draws heavily upon previous work |^ in which we have adopted 
Lagae's method to investigate the nature of full Salpeter solutions. Our prin- 
cipal conclusion was that the only linearly confining potential which yields linear 
Regge trajectories and has variationally stable solutions is a time component Lorentz 
vector. This confirms previous work done for the equal mass case 0. 

In the present work we use the fact that stable solutions exist for the time com- 
ponent vector confinement in order to estimate the range of applicability of the 
reduced Salpeter equation. The desirable properties of the time component vector 
potential in the Salpeter equation does not mean that it should be used as a confine- 
ment potential, since it yields wrong sign of the spin-orbit interaction, disagreeing 
both with QCD and experiment. We also compare solutions to the full equation 
and its reduced version for an equal mixture of scalar and time component vector 
confinement. This type of mixed confining kernel has been recently used in |^ for 
the investigation of weak decays of heavy mesons. The vector confinement stabilizes 
the scalar confining part up to the case of equal mixtures. Phenomenologically, the 
scalar confining part is necessary to reduce the P-wave spin-orbit splitting. For this 
mixed confinement case we also explicitly demonstrate that the reduced Salpeter 
equation is adequate for the investigation of the heavy-light systems, such as D 
and B mesons, as well as for heavy onia. We also examine the extent to which ^ 
corrections to heavy-light systems depend on which wave equation is used. We find 
that the difference is negligible even for D mesons. 

In Section |^ we adapt Lagae's formalism to the reduced Salpeter equation. 
Our numerical results are contained in Section |^ where we compare the full and 
reduced Salpeter solutions for both onia and heavy-light mesons. Our conclusions 
are summarized in Section |. In the Appendix ^ we provide the complete reduced 
Salpeter radial equations for the three Lorentz type kernels, 7°® 7° [time component 
vector], 1 (g) 1 [scalar], and ® 7^^ [full vector]. 



3 



2 Reduced Salpeter equation 



Recently Lagae has proposed an elegant formalism for the reduction of the full 
Salpeter equation to a system of equations involving only radial wave functions. 
One of the nice things about his method is that the transition from full to reduced 
Salpeter equations can be accomplished easily. In this section we briefly sketch the 
main points of this formahsm as adapted to the reduced Salpeter equation. 

We start from the Salpeter equation for a fermion-antifermion system in the CM 
frame of the bound state, 

d^k' rA^(k)70[y(k, k')$(k')]7°A2 (-k) 



$(k) 



M -El- E2 



M{k)-f^[V{k, k')$(k')]7°A^(-k) 



M + Ei + E2 

Here, A!j_'s are the usual energy projection operators, given by 

, E,ik)±H,ik) 
^ 2E,{k) ' 

with Hi being the generalized Dirac Hamiltonians, 

H,{k) = A{k)a-k + B,{k)(3 , 



(1) 



(2) 



(3) 



and Ei{k) = J Ai{ky + i?j(k)^. Again, we'll consider constituent quarks of masses 



rrii, so that 



A(k) 
5.(k) 
E^{k) 



rrii 



+ k2 . 



(4) 
(5) 
(6) 



The formal product of in equation (0) represents the sum of scalar potentials Vi 
and bilinear covariants, 



y(k,k')$(k') -^E^*(k,k')G,$(k')G', , 



(7) 
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where the Gj's are Dirac matrices. 

The reduced Salpeter equation is obtained by dropping the second term from 
(|l|), and this is usually justified for heavy-quark systems on the grounds that 

M - El- E2 

< 1 . (8) 

M + E1 + E2 ^ ^ 

The resulting equation, 

is a standard eigenvalue equation, and it has been used in a number of studies of 
relativistic bound states P, P|, PH . 



Mm = {El + E2)m + J A^(k)7°[y(k, k')$(k')]7°A'-(-k) , (9) 



In order to apply Lagae's formalism to the reduced Salpeter equation, we 
multiply (H) by 7°, and define 

X(k) = $(k)7°, (10) 

Ti = 7°G, , (11) 



so that (|9f) becomes 



Mx = {El + E2)x + E / 7Tr^^^(k - k')A^r,x'r.A2_ , (12) 

where notation / = /(k), /' = /(k') is employed. Vi(k — k') has the Fourier 
transform V{r) in the case of Lorentz vector kernel, and —V{r) in the case of Lorentz 
scalar kernel. 

Using properties of projection operators, it can be easily shown that the full 
Salpeter amplitude satisfies 

^X + X-^ = 0. (13) 

tLi 

For the reduced equation, this constraint breaks into two parts, 

Hix = Eix, (14) 
XH2 = -E2X. (15) 
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Taking these constraints into account, the norm of the reduced Salpeter amphtude 
TT| , ^ can be written as 



\X\ 



(2vr) 



rTr 



and is related to the normahzation of bound states as 

,, ,,o 1 



.(B\B) . 



(16) 



(17) 



Using (|T^ inside of (|T^) one obtains 



(18) 

This equation will be used for obtaining radial equations from the variational prin- 
ciple as outlined in [Q. It is interesting to note that it has the same form for both 
full and reduced Salpeter equations. 

Now, in the case of the full Salpeter equation, one expands the amplitude as 



X = f-o 



Co + Cipi +Afo - a + Afi- 



Picr 



(19) 



with 16 Hermitian matrices whose squares are unity [1, pi, cr, picr) defined in 0. 
Using this decomposition, it is easily seen that constraint ([T3|) can be satisfied 
by expressing the 16 components of x (-^'s and A/''s) in terms of eight functions 
(Li, L2, Ni, N2). The correct form for £'s and A/''s is given in 0. For the reduced 
Salpeter equation, both constraints ( [T^ ) and (|T5|) can be simultaneously satisfied if 
Li = L2 = L and Ni = N2 = N. 



Following 1 13], we obtain the radial equations by expressing L and N in terms 



of spherical harmonics and vector spherical harmonics [14|, so that 



L(k) 
N(k) 



L{k)YjMi'k) , 

iV_(A;)Y_(k) + iVo(A:)Yo(k) + iV+(fc)Y+(k) , 



(20) 
(21) 
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where Y_, Yq, and Y+, stand for Yjj_im, Yjjm, and Yjj+im, respectively. We 
also introduce functions n+ and n_, defined as 







fj, V 










— V fl 







(22) 



with 

Using these definitions inside expressions for the £'s and AT's as given in 0, together 
with properties of spherical and vector spherical harmonics, ([T8|) can be expressed 
in terms of radial wave functions only. Then by taking variations with respect to 
L*{k), NQ{k),n\{k), and n*_{k), as explained in @], one obtains the set of coupled 
equations for the radial wave functions of the reduced Salpeter amplitude. We 
summarize these equations in Appendix ^ for the kernels 7*^ ® 7°, 1 ® 1 and 7^* ® 7^. 

3 Numerical results 

As outlined in Appendix B of [Q, one can solve the system of radial equations by 
expanding the wave functions in terms of a complete set of basis states, which depend 
on a variational parameter /?, and then truncating the expansion to a finite number 
of basis states. In this way, a set of coupled radial equations can be transformed into 
a matrix equation, Tiip = Mi/j. The eigenvalues M of the matrix Ti will depend on 
P, and by looking for the extrema of M{(3), one can find the bound state energies. If 
the calculation is stable, increasing the number of basis states used will decrease the 
dependence of the eigenvalues on (3. Regions of (3 with the same eigenvalues should 
emerge and enlarge. For each of the results discussed in the remainder of this paper 
we have verified that this indeed occurs. 
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3.1 Equal mass case with 7 (g) 7 kernel 

In Figure we compare solutions of reduced and full Salpeter equations for equal 
mass systems with a pure time component vector confinement {V{r) = ar, a = 
0.2 GeV"^). We have varied the quark masses (mi = m2 = m) from to 1 GeV, 
solved both equations for all J = 0, 1, and 2 states (which involves all S, P, and 
most D waves), and plotted the difference between state mass and rest mass of the 
two quarks. As one can see, the difference between the two solutions is noticeable 
only for J^^ = 0~~^ and 0"^"*" states, and then only for very small quark masses. For 
example, for zero quark masses the difference for the O""*" state is about 25 MeV, 
while already for quark masses of 0.3 GeV it is only 6 MeV. On the other hand, 
for the 1 state the difference between the two solutions is about 1 MeV even for 
zero quark masses. Another interesting thing to observe in Figure |I] is that for both 
equations and for zero quark mass we have degeneracy of 0~^ and 0^+, 1 and 
1"'""'", and also 2+"*" and 2 states. This parity degeneracy can be easily explained 
by referring to the radial equations for the full Salpeter equation given in Appendix 
A of In the limit where both masses go to zero, it can be easily seen that 0""*" 
and O"*"^ equations are the same. Similarly, for J > states the four radial equations 
for P = (— l)"^ and G = (—1)"' (involving ni+,n2+,ni_ and n2_) decouple into two 
systems of two equations. The first one (involving ni+ and ^2+) is the same as 
the system describing P = (—1)"'+^ and G = (—1)"^^^ states, while the second one 
(involving and ^2-) is equivalent to the system describing P = (—1)"'+^ and 
G = (—1)"' states (and higher in energy, as can be seen in Figure The m = 
degeneracy is an example of the chiral symmetry of the vector potential and its 
Wigner-Weyl realization through parity doublets. 

In order to see the effects of the short range Coulomb potential, we have per- 
formed a similar analysis with V{r) = ar — ^, using a = 0.2 GeV^ and k = 0.5. 
The results are shown in Figure ^. Again, the difference between full and reduced 
Salpeter solutions is noticeable only for the J^*^ = and 0^+ states. For the 
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0^^ state, the difference is now about 35 MeV for mi = 1112 = m = 0, and about 
10 MeV for mi = m2 = m = 0.3 GeV. For tlie 1 state, tlie difference is only 
about 3 MeV for zero quark masses. 

Finally, in Figure ^ we show the results of the same analysis as above, but this 
time including an additive constant, V{r) = ar + C ~ ^ {a = 0.2 GeV'^, C = 
— 1.0 GeV, and k = 0.5). One might expect that adding a constant to the potential 
would not change the difference between the two equations. However, as one can 
see from the Figure ^, it is not quite like that. Now the solutions to the full Salpeter 
equation for the O""*" and 0''"+ states are considerably lower in energy than the 
solutions to the reduced Salpeter equation. For the 0^+ (1 ) state the difference 
is about 106 MeV (7 MeV) with zero quark masses and about 3 MeV (1 MeV) 
with mi = 7722 = 1-0 GeV. 

The reason for this somewhat unexpected behavior is that a negative constant 
G added to the kernel of the full Salpeter equation lowers the eigenvalues by an 
amount larger than |C|, while for the reduced Salpeter equation it is exactly \G\. For 
example, adding C = — 1.0 GeV to the potential V{r) = ar — ^ with a = 0.2 GeV'^ 
and K = 0.5 the lowest eigenvalue for the 0~"^ state (with zero quark masses) is 
lowered by about 1.072 GeV for the full Salpeter equation with a time component 
vector kernel. This effect is much less noticeable with larger quark masses, and 
higher J states, e.g. for the 1 state with zero quark masses and same a and k as 
before, the lowest eigenvalue was lowered by 1.004 GeV after adding C = — 1.0 GeV. 
We also note that these numerical results were obtained with 25 basis states, so that 
dependence of the results on the variational parameter characterizing the basis states 
was negligible. 

In order to further explore the relationship between the full Salpeter equation 
and the reduced one, we have plotted the radial wave functions for the O""*" case and 
for V{r) = ar + C - f (a = 0.2 GeV^, C = -1.0 GeV, and n = 0.5). Just as a 
reminder, the reduced Salpeter equation for the pseudoscalar case has only one wave 
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function (L), as opposed to two (Li and L2) in the full equation. Also, when the 
reduced Salpeter equation is valid, then Li and L2 are equal. As we can see from 
Figure for very small quark masses (mi = 1112 = 0), the difference between Li and 
L2 is large, and the reduced equation cannot replace the full one. However, with 
nil = = 1.0 GeV (Figure |^), the reduced Salpeter result is much more closer 
to the full one. In these two figures we use a Cornell potential with an additive 
constant (a = 0.2 GeV^, C = -1.0 GeV, and n = 0.5). 

From this analysis, it is clear that the solutions of the reduced Salpeter equation 
are nearly the same as those of the full one for the description of the heavy-heavy (cc 
and bb) mesons, and a very good first approximation even for the ss mesons (with 
s quark mass of about 500 MeV). This justifies the assumption of Gara et. al. |]TU| 
that the reduced Salpeter equation could be used for the description of ss mesons. 

3.2 Heavy-light case with 7^ 7^^ kernel 

A similar analysis can be performed for the "heavy-light" systems. For V{r) = 
ar + C - f , with a = 0.2 GeV^, C = -1.0 GeV, and k = 0.5, we fixed the "light" 
quark mass at mi = 0, varied the "heavy" antiquark mass m2 from to 1 GeV, 
and solved both equations for all J = 0,1 and 2 states. The results are shown 
in Figure ^ The degeneracy of states with the same J and different parity can 
be again explained easily by looking into the radial equations for the full Salpeter 
equation given in Appendix A of In the limit where mi — > 0, 0i — * and 
(p —>■ 9, which makes equivalent the two sets of equations for different parities. As 
far as the difference between the full and reduced Salpeter equations are concerned, 
it is again important only for J = states. For example, for 0~ and O"*" (1~ and 
1+) it is only about 7 MeV (1 MeV) at m2 = 1.0 GeV. Figure shows that for 
mi = 0, m2 = 1 GeV, L is already a very good approximation to Li and L2. One 
also has to remember that with such a large negative constant the c quark mass 
must be considerably larger than 1.0 GeV in order to describe D mesons. Given all 
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this, we conclude that the reduced Salpeter equation is an excellent approximation 
to the full one for the description of D and B mesons. 

Although the time component vector interaction has many nice properties, it 
is flawed as a realistic quark confinement interaction. As pointed out earlier, it 
predicts "parity doubling" of meson states in the limit of zero quark mass. For large 
quark masses this difficulty appears as the "wrong sign" spin-orbit interaction which 
conflicts both with experiment and QCD. 

3.3 Mixed confinement potentials 

As already mentioned, recently a half-half mixture of the time component vector 
and scalar confinement has been proposed in , together with a one gluon exchange 
kernel, for the investigation of weak decays of B and D mesons. In order to compare 
the full Salpeter equation with its reduced version in this type of model, we adopt 
the mixed confining kernel, 

^[7°®7° + l® l]K(r) , (24) 

with 

V^{r) = ar + C , (25) 
and for the short range potential we simply take 

[7°®7°]V;(r) , (26) 

where 

V,{r) = . (27) 

A confinement mixture of this type has been shown to have a stable variational 
solution p. For the parameters of the potential we choose a = 0.2 GeV^, C = 
—1.0 GeV, and n = 0.5. Computation of the equal mass case is shown in Figure 
|] (for the 0~"^ and 0^+ states). As one can see, the differences between full and 
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reduced Salpeter solutions are only slightly different than in the case with a pure 
time component vector kernel. For the 0"^ state and rrii = m2 = m = 1.0 GeV 
the difference between the two equations is about 7 MeV. The heavy-light case 
calculation (for the same potential parameters) is shown in Figure ^ The difference 
between the two solutions for the 0" state, and for mi = and m2 = 1.0 GeV, is 
about 9 MeV. Therefore, we again conclude that the reduced Salpeter equation is 
as good as the full Salpeter equation for the description of the cc and bb mesons, a 
very good first approximation even for the ss mesons, and would serve as well as 
the full Salpeter equation for the description of the heavy-light systems, such as D 
and B mesons. 

Of course, these results are dependent on parameters of the particular model. 
However, in our analysis we have used values for a and k, that are typical in the 
hadron spectroscopy, and constant C that is perhaps slightly larger than usual. We 
have also restricted ourselves to constituent masses that are smaller than the usually 
assumed c quark mass. Therefore, we feel that our main conclusions would not be 
drastically altered if a different set of realistic parameters was used. 

In order to illustrate this, we have chosen parameters of the potential to be as 
close as possible to the ones used in (as given in their Table 1), i.e. mi = 0.2 GeV, 
ma = 1.738 GeV, a = 0.335 GeV^, C = -1.027 GeV, and k = 0.521 (which 
corresponds to agat = 0.391 in and solved both equations for the 0~ and 0"*" 
states, with the kernel described by ( p4| - ^ ). The differences between ground state 
energies were 5 MeV and MeV, respectively, despite the large value of a. For the 
0~ state, where the difference between the two solutions should be most obvious, 
we have plotted the radial wave functions in Figure |TD|. As one can see, the reduced 
wave function is a very good approximation for the full wave functions. 

For the sake of simplicity, in the previous calculations we have used a short range 



potential with a fixed coupling constant, for which Murota has shown most of 
the Salpeter amplitudes are divergent as r ^ 0. If one uses a running coupling 
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constant, this divergence is less pronounced, but still present. That is precisely the 
reason why the short range potential used in was regularized. In order to show 
the effects of regularization, instead of (|27|) we now take as in 

{4 a(r) 
-3—' ^-^V (28) 
agT^ + bg , r < ro 

The constants ag and bg are determined by the condition that Vg{r) and its derivative 
are continuous functions. The running coupling constant is parametrized exactly as 
in 0, with their value of ro = 0.507 GeV^^, and the saturation value for the 
coupling constant asat = 0.391. The string tension and constant were again a = 
0.335 GeV^ and C = —1.027 GeV, and quark masses were mi = 0.2 GeV and 
1712 = 1-738 GeV, as for the previous calculation. Using these parameters, we have 
again solved both equations for the 0~ and 0"*" states. The differences between ground 
state energies were 3 MeV and MeV, respectively, showing that a regularized short 
range potential reduces the differences between the reduced and the full Salpeter 
equation. For the 0^ state, we have again plotted the radial wave functions in 
Figure |TT]. As one can see, all wave functions are now finite at the origin, and the 
reduced Salpeter wave function is an even better approximation to the full ones than 
it was before. 

We can also use this model to estimate the accuracy of — recoil corrections to 
the heavy-light limit. In Figure |T2| we show the difference between the 0~ ground 
states for a finite and an infinite heavy mass (rri2) with a massless light quark (mi) 
in both cases. We see that these "recoil" corrections are quite important even for 
the b quark mesons where correction is nearly 40 MeV (at ^ 0.2 GeV^^). On 
the other hand, the difference between full and reduced Salpeter solutions is small. 
For a charmed meson (— ~ 0.66 GeV~^) the difference is about 3.5 MeV, while 
for a meson with a b quark it is about 0.2 MeV . 

In the mixed confinement ( plf ) was used in part because the full Salpeter 
equation does not have stable solutions unless the scalar confinement part is equal 
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or less than the time component vector part. We should note that the pure scalar 
confinement could have been used with the reduced Salpeter equation. 

4 Conclusions 

The reduced Salpeter equation, also known as the no-pair equation, has long been 
used in dynamical models of mesons. It has also been long appreciated that it is 
an approximation to the full Salpeter equation and that the discarded portion only 
vanishes if at least one of the constituent masses is infinite. The reduced equation 
has nevertheless been used because it has the standard hermitian form. 

In this paper we have examined the conditions under which the reduced equation 
can be employed without significant loss in accuracy. The critical factors turn out 
to be constituent mass, state, and the nature of the interaction. If the total 
quark mass exceeds about 1.0 GeV very little difference is found between the full 
and reduced Salpeter solutions. Also, with the exception of the 0~ and 0"*" states 
very little difference is found even at zero quark mass. Finally, even for 0^ states and 
vanishing quark mass the differences between full and reduced Salpeter solutions are 
small if there is no large constant in the coordinate space confining potential. 

There remain a number of hadronic states with light quark masses in which the 
full Salpeter equation must be used. Differences up to 100 MeV were found between 
pseudoscalar masses at zero quark mass for the two equations. 

In our comparison between the full and reduced Salpeter solutions we have con- 
sidered both energies and wave functions. As was the case with the energy eigen- 
values, we see large differences between the 0~ full and reduced wave functions for 
zero quark mass (see Fig. The differences are largest at the origin, r = 0. As 
observed in subsequent figures, increasing the quark mass and considering higher 
states causes the reduced Salpeter wave functions to become more similar to the 
fuU ones. The difference between the two solutions is always most noticeable at the 
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origin. 

We have primarily considered the time component vector kernel, since its solu- 
tions with the full Salpeter equation are variationally stable and yield normal linear 
Regge trajectories in the case of linear confinement 0. Although the solutions with 
a time component vector potential have many desirable properties, a quark confine- 
ment of this type has a spin-orbit interaction of the wrong sign. The addition of up 
to equal parts Lorentz scalar confinement has been advocated recently |^ for the 
study of weak decays of heavy-light mesons. The variational stability is retained in 
this case and the reduced Salpeter equation is shown to be accurate under similar 
conditions as in the pure time component vector case. The reduced equation has 
the additional advantage of variational stability with pure scalar confinement. 
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APPENDIX 



A Radial equations 

In this appendix we give the final form of the radial equations for the reduced 
Salpeter equation for the kernels 7° 7°, 1 ® 1 and 7^ CS) 7^. These equations 
represent a general case with a quark of mass mi and an anti-quark of mass m2- 
However, one has to keep in mind that for J = two wave functions vanish, i.e. we 
have A^o = 0, and n+ = 0. 

As in p, 1^ we have used notation 

Sff, = sin , = cos (j) , (29) 
Sg = sin 6* , Cg = cos 9 , (30) 



with angles (p and 9 defined as 

01 + 02 n 02 - 01 



(31) 



while 0j's are defined through 



A- B 

cos 0j = — ^ , sin 0j = . (32) 



Ai, Bi and Ei are defined in 

In the equal mass case the equations given below somewhat simplify, since one 
has El = E2, = 01 = 02, and 6' = 0, so that Sg = and Cg = 1. Also, since 
charge conjugation is a good quantum number in the equal mass case, the two P = 
(— 1)'^^^ state equations decouple into two separate equations, one corresponding to 
C = (—1)"' (involving L), and the other corresponding to C = (—1)"^^^ (involving 
No). 

The heavy-light limit {m,2 00) is obtained by setting £"2 — ^ m,2, 02 — > f , so 
that Sg Ca) and Cg ^ Sa,. As expected, in the heavy-light limit equations for the 
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7° ® 7° and 7^ 7/^ kernels are the same. Also, in this case spin of the heavy quark 
decouples from the spin of the light quark, so that total angular momentum j of the 
light quark becomes a good quantum number. Inverting (^21) , 



= lifi^ — un_ , (33) 

= ijn+ + fm_ , (34) 

and also putting 

= i/L - fiNo , (35) 

L_ = //L + uNo , (36) 

from the heavy-light limit equations in terms of n+, n_, Nq and L, one can obtain 
decoupled equations in terms of A^+, A^_, and L_, describing heavy-light states 
with quantum number j. There will always be a pair of degenerate states, described 
with and L+ ( J = L + 1 and J = L, for the state with j = L + and and 
L_ {J = L — 1 and J = L, for the state with with j = L — 

For any mixture of different kernels, only the kernel parts of the radial equations 
should be added. The kinetic energy terms are always the same. In the 1 ® 1 case, 
we have introduced an additional minus sign in the kernel, so that V{r) has the same 
form for all three cases considered, e.g. for the Cornell potential V{r) = ar — ^. 

A.l 7^ 7^ kernel 

States with parity P = (—1)'^+^: 

ML = \E, + E^]L + -j^ iCeVjC'eV + S^VjS'^L' 



+ ^llySe{V.J^l - Vj+i)C'^N', + /i//C^(y7_i - Vj+i)S'sK] , (37) 

1 k''^dk' 

MNo = [E^ + E2]No + - -—^[CeVjC'sN;, + S^VjS'^N;, 

2 Jo (iTTj^ ^ 
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+ /xz/C<^(yj_i - Vj+i)S',L' + iipSe{Vj.r - Vj+^)C'^L'] . 
States with parity P — (—1)"^: 

Mn+ = [E^ + E2]n+ + -j^ -^[C^VjC'^n'^ + SeVjS'eU^ 

+ iivCeiVj.^ - Vj+i)S'^n'^ + //^^^^(yj-i - Vj+i)C',n'^] . 

A. 2 10 1 kernel 

States with parity P = 

ML = [E^ + E,]L + -J^ ^[CeVjC'eL' + S^VjS'^L' 

- Seifi^Vj., + p^Vj+i)S',L' - C<^(/iVj_i + z/Vj+i)C;l' 

1 POO y'^rlh' 

MNo = [E, + E2]No + -J^ j^[CeVjC'gN;, + S^VjS'^N;, 

- fiuC^iVj^, - Vj+ML'] - fiuSg{Vj_, - Vj+i)C'^L'] . 
States with parity P = (—1)'^: 

Mn+ = [E, + E,]n+ + -J^ ——[-C^VjCy_,-SeVjS'sn'^ 
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+ piuS^iVj^i - Vj+,)C'gn'_ + f^uCeiVj.i - Vj+i)S'^n'_] , (40) 
= [Ei + E2]n. + -j^ -j^[-C^VjC'^n'_-SeVjS'en'_ 



A. 3 kernel 

States with parity P = 

ML = [E, + E2]L + -^[2CbVjC',L' - S^VjS'^L' 

+ Seifi'^Vj.i + p'Vj+i)S'gL' + fxuSeiVj.i - Vj+,)C'^N'^ , (41) 

+ C^{u^Vj.^ + ^Vj+i)C;iV^ + iivC^iVj.r - Vj+,)S'sL' 
States with parity P — (—1)"^: 



roc h n h 



Mn^ = [E, + E,W + j^ 

+ C,(z.Vj_i + //Vj+i)C^nV + ^JivCe{VJ-l - Vj+,)S'^n'_] , (42) 

POO y^fjh' 
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FIGURES 



Figure 1: Equal mass comparison of the reduced (dashed Unes) and full Salpeter 
(full lines) equations for the time component vector kernel with V{r) = ar {a = 
0.2 GeV'^). The energies of all states with J equal to 0, 1, and 2 are shown as a 
function of the quark mass. We have used 15 basis states. 

Figure 2: Equal mass comparison of solutions to the reduced (dashed lines) and full 
Salpeter (full lines) equations for the time component vector kernel with V{r) = 
ar — ^ (a = 0.2 GeV'^, k = 0.5). The energies of all states with J equal to 0, 1, and 
2 are shown as a function of the quark mass. We have used 15 basis states. 

Figure 3: Equal mass comparison of the reduced (dashed lines) and full Salpeter 
(full lines) equations for the time component vector kernel with V{r) = ar + C — ^ 
{a = 0.2 GeV"^, C = —1.0 GeV, k = 0.5). The energies of all states with J equal to 
0, 1, and 2 are shown as a function of the quark mass. We have used 15 basis states. 

Figure 4: Pseudoscalar {J^'~'' = 0^"*") radial wave functions in coordinate space for 
the reduced (L, dashed line) and full Salpeter equations (Li, lower full line, and 
L2, upper full line), with time component vector kernel and V{r) — ar + C — ^ 
(a = 0.2 GeV'^, C — —1.0 GeV, k — 0.5). The quark masses were rui — m2 — 0, 
and the calculation was done with 25 basis states. 

Figure 5: Pseudoscalar (J^*-^ = O""*") radial wave functions in coordinate space for 
the reduced (L, dashed line) and full Salpeter equations (Li, lower full line, and 
L2, upper full line), with time component vector kernel and V{r) — ar + C — ^ 
{a = 0.2 GeV^, C = -1.0 GeV, k = 0.5). The quark masses mi = ma = 1 GeV , 
and the calculation was done with 25 basis states. 
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Figure 6: Comparison for heavy-light mesons of the reduced (dashed hues) and full 
Salpeter (full lines) solutions for the time component vector kernel with V{r) = 
ar + C - ^ {a = 0.2 GeV^, C = -1.0 GeV, k = 0.5). The lighter quark mass 
was fixed at mi = 0, and we show the light degree of freedom energy M — m2 as a 
function of m2 for the lowest angular momentum states . We have used 15 basis 
states. 

Figure 7: Heavy- light pseudoscalar (J^ = 0~) radial wave functions in coordinate 
space for the reduced (L, dashed line) and full Salpeter equations (Li, lower full line, 
and L2, upper full line), with time component vector kernel and V{r) = ar + C — f 
(a = 0.2 GeV"^, C = —1.0 GeV, k = 0.5). The quark masses were mi = and 
7722 = 1 GeV. The calculation was done with 25 basis states. 

Figure 8: Equal mass comparison of the reduced (dashed lines) and full Salpeter (full 
lines) ground state O""*" and O"*"^ energies. An equal mixture of the time component 
vector and scalar confinement (V^(r) = ar + G), together with time component 
vector short range potential {Vg{r) = — ^) was used. The potential parameters were 
a = 0.2 GeV^, C = —1.0 GeV, and k = 0.5. Comparison with Figure § shows the 
breaking of the parity degeneracy at m = 0. 15 basis states was used for calculation. 

Figure 9: Heavy-light mixed confinement comparison of the reduced (dashed lines) 
and full Salpeter (full lines) ground state 0~ and 0"^ energies. An equal mixture of 
the time component vector and scalar confinement {Vc{r) = ar + C) , together with 
time component vector short range potential (V^(r) = — ^) was used. The potential 
parameters were a = 0.2 GeV^, G = —1.0 GeV, and k = 0.5, while the lighter 
constituent mass was fixed at mi = 0. By comparing with Figure ^we observe the 
lifting of the parity degeneracy present in a pure time component vector interaction 
at m = 0. 15 basis states was used in the calculation. 
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Figure 10: Pseudoscalar {J^ = 0^) radial wave functions in coordinate space for 
the reduced (L, dashed hne) and full Salpeter equations (Li, lower full line, and 
L2, upper full line), with a half-half mixture of the time component vector and 
scalar confinement (V^(r) = ar + C), together with time component vector short 
range potential {Vg{r) = — The potential parameters were a = 0.335 GeV"^, 
C = —1.027 GeV, and k = 0.521, while the quark masses were nii = 0.2 and 
m2 = 1.738 GeV. The calculation was done with 25 basis states, and represents a 
model of but with a singular short range potential. 

Figure 11: Pseudoscalar {J^ = 0^) radial wave functions in coordinate space for 
the reduced (L, dashed line) and full Salpeter equations (Li, lower full line, and L2, 
upper full line), with a half-half mixture of the time component vector and scalar 
confinement (V^(r) = ar + G), together with the regularized time component vector 
short range potential (as defined in (|2^) in the text). The potential parameters 
were a = 0.335 GeV^, G = -1.027 GeV, a^at = 0.391, and tq = 0.507 GeV-\ while 
quark masses were mi = 0.2 and 1712 = 1.738 GeV. The calculation was done with 
25 basis states, and represents a model of |0, including a regularized short range 
potential. 

Figure 12: ^ corrections to the heavy-light 0~ ground state energy as a function of 
^ using the model of 0, with light quark mass mi = 0, and potential parameters 
the same as before. The correction ranges from about 40 MeV for the B meson 
to about 100 MeV for the D meson. The difference between the full and reduced 
solutions is 0.2 MeV and 3.5 MeV for the B and D mesons respectively. 
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